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Abstract: The Sound attenuation and dispersion in saturated gas-vapor-droplet mixtures with 
evaporation has been investigated theoretically. The theory is based on an extension of the work 
of Davidson (1975) to accommodate the effects of transpiration on the linear particle relaxation 
processes of mass, momentum and energy transfer. It is shown that the inclusion of transpiration 
in the presence of mass transfer improves the agreement between the theory and the experimental 
data of Cole and Dobbins (1971) for sound attenuation in air- water fogs at low droplet mass 
concentrations. The results suggest that transpiration has an appreciable effect on both sound 
absorption and dispersion for both low and high droplet mass concentrations. 
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= vapor mass loading (ratio of mass of vapor to mass of gaseous mixture per unit 
volume of undisturbed mixture-droplet system, 

- droplet diameter 

= mass diffusivity of the gaseous mixture 
= frequency of sound 

= drag force between the mixture and a single droplet, — - m)^ 

8 

= convective heat transfer coefficient between the mixture and a single droplet 
= latent heat of vaporization, )(rp - To ) 

= convective mass transfer coefficient between the mixture and a single droplet 
= specific heat ratio of gas, 

= specific heat ratio of particles, 

= specific heat ratio of particles, 

= mixture thermal conductivity 
= complex wave number, + ik^ 

= latent heat parameter (dimensionless), /i^o / R^Tq 
= Lewis number, k / pc^D 

= mass of one droplet 

= diffusional mass transfer between the mixture and a single droplet, wdpDSh{p^ ~ P,) 

= mean number of droplets per unit volume of mixture 
= droplet Nusselt number, hdp / k^^ 

= mixture pressure 
= Prandtl number, Cp^p^ / k,^ 

= heat transfer between the mixture and a single drop, Tid pkNui^p -T) 

= ideal gas constant of mixture (instantaneous), Rq (} + R') 

= droplet Reynolds number, /Oq|m^ - u^dp / 

= droplet Sherwood number, h^dp / Dq 

= Schmidt number, = Le Pr 

= Spalding transport parameter 
= mixture temperature 
= mixture velocity in the x - direction 
= coordinate distance 


Greek Symbols 

a = sound (energy) absorption coefficient (dimensionless) 

P = sound dispersion coefficient, (cq / c^)- 1 

Y - mixture isentropic exponent (ratio of specific heats), 

= mixture dynamic viscosity 

p = mixture density (mass per unit volume of gaseous mixture) 

= droplet mass per unit volume of mixture-droplet system 
- particle mass transfer relaxation time, Pp I llD^p^ = (3 / 2)Sctj 
= particle dynamic (momentum) relaxation time, d^pPp !\%p 
r, = particle thermal (heat transfer) relaxation time, Pr Cppd^ Pp li^lp^Cp^ )=(3/2)Prr^ 

Q) = circular frequency of sound 

Subscripts 

0 = property in undisturbed medium 

1 = transpiration 

g = gas (nonvolatile component, air) 

p = droplet, particle 

s = saturated vapor 

V = vapor 

1 . Introduction 

The propagation of sound (attenuation and dispersion) in a gas containing suspended droplets 
is of considerable interest in many important technical applications such as liquid propellant 
combustion chamber instability [1,2], and noise mitigation by water injection [3-7]. When a 
system comprised of a suspension of small spherical liquid droplets in equilibrium with a 
gaseous mixture is perturbed by a plane acoustic wave, fluctuations arise in the velocity, 
pressure, temperature and density of the gaseous mixture. Oswatitsch [8] treated sound 
propagation in a gas containing volatile particles exchanging heat (thus allowing 
evaporation/condensation), but excluded the particle viscous drag. In their classical work. Marble 
[9] and Marble and Wooten [10] treated sound attenuation in a saturated vapor-gas-droplet 
mixture (small spherical droplets initially in thermal equilibrium with the gaseous mixture) 



considering mass diffusion, momentum and heat transfer. They neglected differences between 
the gas and vapor properties, and assumed that the gaseous phase is both thermally and 
calorically perfect when mass transfer is present [11,12]. They also incorrectly dropped a term in 
the energy equation (for the mixture enthalpy) during linearization [11]. 

Cole and Dobbins [13] theoretically studied the attenuation of a plane acoustic wave with 
specific reference to air-water fogs (low droplet mass fraction), accounting for the difference 
between the gas constants for the vapor and the noncondensible gas, and incorporating the 
correct form of the energy equation for the mixture enthalpy. But they incorrectly defined the 
mixture enthalpy, and neglected the fluctuations of the gas constant of the inert gas-vapor 
mixture (retained the assumption of calorically perfect gas mixture). Their theoretical predictions 
deviate (overpredict) by about 35 percent from their measurements [14] for sound attenuation in 
atmospheric fogs (Fig. 1), largely due to the failure to account for the fluctuation in the mixture 
gas constant in the presence of mass transport. 

Davidson[ll] accounted for the individual gas constants for the vapor and the 
noncondensible gas (air), dispensed with the assumption that the gaseous mixture is calorically 
perfect, and incorporated the correct form of the linearized energy equation. With these 
improvements, better agreement between the theory and the experimental data of Cole and 

Dobbins [14] is achieved. However the agreement still fell short as evident from Fig. 1. 

\ 

Gubaidullin and Ivandaev [12] took into account the transient and nonequilibrium effects of 
interphase mass, momentum, and energy (heat) transfer on the propagation of acoustic waves in 
gas-droplet mixtures. They concluded that the nonequilibrium phase transition theory has only a 
slight (not appreciable) improvement as the nonequilibrium curves shift the attenuation curves to 
the left without changing the maximum absorption coefficient appreciably. Nonequilibrium 
considerations with regard to the absorption and dispersion of sound were also considered by 


Gumerov et al. [15]. 



From the theory of convective heat and mass transfer, it is well known that transpiration 
tends to reduce droplet drag, and heat and mass transfer (Bird et al. [16]). Dupays [17] 
considered the transpiration effect on sound propagation in a droplet-gas mixture considering the 
general case that the droplet is in the thermal nonequilibrium with the gaseous mixture. The 
conservation equations of the gas-droplet mixture are essentially those of Davidson [11] except 
for a minor difference in the energy equation. However, their theory (either equilibrium or 
nonequilibrium) was not validated with any experimental data concerning the role of 
transpiration on sound attenuation. 

The present work endeavours to investigate the role of transpiration on the attenuation and 
dispersion of sound in evaporating gas-droplet mixtures initially under thermodynamic 
equilibrium. The approach is based on an extension of the theory of Davidson [11] to 
accommodate transpiration on the linear relaxation processes. The results of the transpiration 
theory will be validated by a comparison of the predictions with the only available experimental 
data of Cole and Dobbins [14] for sound attenuation in atmospheric fogs. The effect of 
transpiration in mixtures characterized by relatively high droplet mass concentration (neglecting 
particle-particle interaction effects) will then be investigated with regard to sound attenuation 
and dispersion. 

2. Extension of Davidson's Theory for transpiration 

Considering plane wave solutions of the form 

<t>! <t>^- exp[/(/f<m ! CQ-coif[ . ( i ) 

Davidson [11] arrived at an expression for the dimensionless complex wave number for plane 
wave propagation as follows: 

K = k^+ik2 = f(T^,Tt,Tc,C„,C^,Hp,H^,Hg,?r,Le,r,rv,L,R^/Rg) ( 2 ) 
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In the preceding equation, the quantities denote respectively the particle relaxation 

times for momentum, heat and mass transfer: 
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Physically the quantities , r, , represent a measure of the inability of the particles to follow 
(respond to) the fluctuations in the fluid motion, heat and mass diffusion. The quantity Cp 
denotes the mixture specific heat at constant pressure, pthe mixture density, //the mixture 
viscosity, k the mixture thermal conductivity, D the mass diffusivity, and d the diameter. The 
subscript 0 denotes the undisturbed state, and p the particle. 



The Prandtl number Pr ^ Lewis number Zeand the Schmidt number Sc of gas are defined 
by 

Fr = CpQ^Q/kQ, Le^kl P qc'pqDq, Sc = Vq / Dq = Le?r (b) 

where v = ////? designates the mixture kinematic viscosity. The quantities denote 

respectively the droplet loading and vapor mass loading (saturation condition) defined by 

= p„o / Po = I p^, C, = / Po = C,, (To ) (4c) 

where T denotes the temperature, and the subscript s denotes the saturation state, and v for vapor. 
The remaining parameters are defined by 

H p — c pp / CpQ, Hy — c py / CpQ, L — I RyT(j (4d) 

where c^^is the specific heat of droplet at constant pressure, /z^the latent heat of vaporization, 

and i?the gas constant. The energy attenuation coefficient or and the propagation speed care 
related to the complex wave number K by the relation 

a = 2k2, c/cq = 1/ (5) 

The present analysis extends the work of Davidson [11] to accommodate the effect of 
transpiration (so-called blowing or mass injection into the boundary layer surrounding the 
droplet) on the mass, momentum and heat transfer between the droplet and the surrounding gas. 
Without any loss of generality the propagation of sound in the presence of transpiration may be 
expressed by a relation similar to Eq. (2) as follows: 

K = K{t,^Tp,tJ ( 6 ) 

where the quantities correspond to the relaxation times in the presence of 

transpiration. They are related to the relaxation times in the absence of transpiration by 


the relations 


- ^t\=¥t^t’ ^c\=¥c'^c 


(7a) 


where iyj=CQ^IC[^, y/^=Nu^!Nu, y/^=Sh^lSh , (7b) 

with Cj 3 , , Nu^ , S\ stand respectively for the drag coefficient, Nusselt number and Sherwood 
number in the presence of transpiration. In the absence of transpiration, we have 

Co=24/Rep, Nu = 2, Sh = 2 (8) 

corresponding to the Stokesian drag law, pure heat conduction limit, and pure mass diffusion 
limit respectively. Note that Eq. (8) is valid for low droplet Reynolds number (Re^ «1; 

practically for Re ^ less than about 0.1). Similar modification to the particle relaxation time has 

been incorporated by the author [18, 19] to accommodate the nonlinear drag law and heat 
transfer. 

From the theory of convective heat, mass and momentum heat transfer (Kays and Crawford 
[20]), the net effect of blowing (mass injection) into the boundary layer is to reduce droplet drag, 
heat and mass transfer. The factor y/^ in the presence of blowing (transpiration) is expressed by 

(Eisenkalm et al. [21]) 

Wd{s) = CDxlCo=\l{\ + S) (8a) 

where the quantity 5 is defined by (Kays and Crawford [20]) 

S = {pvs-Pv)^{p~Pvs) (8b) 

and is usually referred to as Spalding's transport number. The following relations for hold 
in the presence of blowing based on a Couette flow approximation (Kays and Crawford [20]) 

y/, {S) = Nu^ t Nu = ln(l + 5)/ S, y/^ (S') = Sh, / Sh = ln(l + S)l S (9) 

The appropriate value of S typical for gas-vapor-droplet mixtures will be determined by a 
comparison of the theory with the experimental data of Cole and Dobbins [14]. 



From physical considerations, transpiration is expected to alter the drag, heat and mass 
transfer which enter into the droplet relaxation times. Its effect on other parameters such as 
mixture composition (less evaporation rates due to transpiration) is believed to be of secondary 
importance to the present problem. 

3. Results and comparison 

3.1 Comparisons with experimental data of Cole and Dobbins [14] 

Figure 2a. presents a comparison of the theory with the experimental data of Cole and 
Dobbins [14] showing the effect of transpiration parameter 5 on the absorption coefficient. The 
data are obtained in a Wilson cloud chamber in which a cloud of nearly uniform size droplets is 
generated. Droplet size and concentration are determined optically with the aid of Mie scattering 

j 

theory for single-size droplets. The size and concentration of droplets are varied over a wide 
range of eontrolling the addition of condensation nuclei [14]. The average drop diameter ranged 
from 1.8 //wto 10/iw (which is about ten times larger than that in atmospheric clouds), the drop 
mass fraction is of the order of 0.01, and the mixture temperature ranged from 271 K to 281 K. 
The frequency of the disturbance is 80 Hz. The range of average drop diameter and droplet 
coneentration result in the value of cot, /C„, ranging from 0.4 to 16 (see Fig. 1). 

The scatter (uncertainty or noise) in the attenuation data is about 1 5 percent in view of the 
difficulties associated with the optical measurement system. To the author's knowledge, these are 
the only quantitative measurements concerning the influence of phase transformation on sound 
absorption in air-water droplet system (see also Baudoin et al. [22]). This error margin is smaller 
than the 35 percent mismatch observed between the theory and the data of Cole and Dobbins 
[14]. 

For the present comparison, the test conditions considered correspond to those indicated by 
Davidson [1 1] in their comparisons corresponding to an average gas temperature of Tq = 276 K. 


=0.00763, Q =0.00807, Pr = 0.715, Le = 0.839 
r = lA, rv=l-32, = 1.61, //g = 1,//, =1.9, //^= 4.2, L = 19.7 

These parameters describe mean fog conditions in the experiment of Cole and Dobbins [14], see 

Davidson [11]. Data points are actually measured over fogs with varying values of C^,C^and 

L . Only data on attenuation were measured, and no data for sound dispersion were obtained. 

Referring to Fig. 2a, the data shown cover a range of coz^ / from 0.4 to 16. The curves for 

5 = 0 correspond to the case of no transpiration. It is seen that with an increase in the value of S 

the peak attenuation (and attenuation in general) decreases, and the peak frequency increases. 

The maximum change in the attenuation is noticed near low frequencies characterized by 

coTf I » 1 . The attenuation at this frequency is primarily due to phase transition (involving 

coupled mass and energy transfer), with the viscous and heat conduction effects being 

unimportant [9]. The effect of S on attenuation becomes relatively small at high frequencies, 

characterized by coT( t in excess of about 4. At 5 = 0.2 , the agreement between the theory 

and the data appears to be most satisfactory relative to the theory of Davidson [1 1] in the absence 
of transpiration (S=0). Thus a value of 5 = 0.2 may be regarded as representative of saturated 
gas-droplet-vapor system. 

Figure 2b. shows a plot of the spectral attenuation as a function of <yr^to reflect the effect of 
transpiration on attenuation due to viscous drag and heat conduction. It is seen that for the 
particular values of and C^, considered, the peak attenuation due to phase transition (coupled 

mass and heat exchange) is considerably high (an order of magnitude larger) relative to that 
arising from viscous and heat conduction effects. On account of the low values of , the two 

peaks are well separated. The values of (ox^ at which these two distinct peaks prevail are clearly 
indicated: the peak due to phase change is seen at (oz^ = , and that due to viscous and heat 



conduction effects at cot^ w 1 (Marble and Wooten [10]). For , which is usually the case, 

the two peaks due to momentum and heat transfer nearly coincide. 

Figure 3. illustrates the dependence of the dispersion coefficient /3 = [(cq/c)^ -l]on the 
transport parameter S . The predictions suggest that the effect of S on the dispersion coefficient 
at low frequencies is appreciable. For cot^ / exceeding about unity, the effect of transpiration 
is not appreciable. It is remarkable that the strongest effect of S on the dispersion is evident at 
the low frequency limit (rur, ->0), whereas the strongest effect of Son attenuation is 

manifested at coTi / = 1 . Unfortunately no experimental data for dispersion are available to 

validate the present theory. 

3.2 Predictions for high droplet concentration 

Calculations are performed for sound propagation in air-water droplet mixtures at relatively 
high droplet mass concentration (C„ ranging from 0.1 to 3) to reflect a range of circumstances in 

typical practical applications involving air- water system. Since the density of the droplet (water) 
is much greater than the surrounding gaseous mixture, the volume fraction of the liquid may be 
neglected compared with that of the gaseous species (Cole and Dobbins [14]; Temkin and 
Dobbins [23]). Thus interaction between neighbouring particles may be neglected even for 
relatively large values of the droplet mass fraction. A vapor mass loading of =0.1 

corresponding to ro=327 K (at a pressure of 0.1 MPa) is considered here. The relevant 
parameters are as follows: 

Pr = 0.7, Le = 0.85, y = 1.4,^ = l.32,R^/Rg -\.6\,Hg = l,H^= 2, Hp = 4.15,L = 15.7 
Figure. 4a shows the predicted attenuation with and without transpiration for =0.1, 0.3, 
1.0 and 3.0. Two distinct peaks in attenuation are noted in the absorption spectrum for =0.1: 
one corresponds to phase transition, and the other due to viscous drag and heat conduction 


effects. For = 0.1 , the attenuation peak due to phase transition exceeds that due to viscous 
drag and heat conduction effects; but the deviation between the two peaks is smaller than that 
corresponding to = 0.01 (Fig. 2b). For values of =0.3 and higher, only one distinct peak 
is noticed, suggesting the existence of interaction between mass, momentum and heat transfer. 
The attenuation peaks increase with an increase in C^, as is to be expected. The peak 
frequencies associated with the viscous drag for both = 0.3, 1.0 and 3.0 are stationed in the 
neighborhood of rur, « 1 , but seem to increase with an increase in . 

With regard to the role of transpiration, the results indicate that transpiration reduces the 
attenuation peaks due to phase exchange and due to momentum and heat conduction. With 
increasing , the reduction of peak attenuation in the presence of transpiration diminishes. At 

= 3 , the peak absorption in the presence of transpiration is slightly enhanced. It is observed 

that the peak frequency increases in the presence of transpiration for all values of droplet 
concentrations considered. Thus the attenuation spectrum is shifted to the right as a consequence 
of transpiration. 

Figure 4b. demonstrates the role of transpiration on the corresponding dispersion coefficient 
for various values of . The results suggest that the dispersion coefficient increases with an 

I 

increase in the droplet concentration , as is to be expected. For frequencies characterized by 
rur,less than about 0.3, the dispersion is relatively insensitive to transpiration. Transpiration is 
seen to have an appreciable effect on the dispersion of sound for uir, »(9(1). Transpiration 
effects are seen to be very small at high frequencies (u)r,>10). The stepwise variation in /?for 
the case of =0.1 corresponds to the separate effects of mass transfer and viscous drag effects 


on the dispersion. 


4. Discussion 


The present theory accounts for transpiration considering only monodispersion of the 
droplets (uniform droplet size). Another possible important mechanism affecting sound 
attenuation is the polydispersion of the droplets (droplet size distribution). In recent times 
Gubaidullin and Nigmatulin [24] and Baudoin et al. [22] investigated polydispersion effects on 
sound propagation, and show that it is important to take into account polydispersion. The results 
of Baudoin et al. [22] suggest that the assumption of monodispersion tends to overestimate the 
absorption coefficient above the peak value (as observed in Fig. 2). The assumption of a mean 
radius tends to overestimate the attenuation, because the influence of large diameters is 
underestimated and the total surface area of droplets is overestimated. Both of the preceding 
articles recommend an average diameter (ij,, based on a mean surface defined as the mean 

volume divided by the mean radius, is satisfactory for both attenuation and dispersion. 

Thus strictly speaking, in polydispersed flow both .transpiration and polydispersion effects 
need to be considered. In the present comparisons, the data pertain to uniform droplet size, and 
thus polydispersion effects are absent. Thus the exclusion of polydispersion mechanism in the 
present development does not invalidate the theory proposed. It is relatively simple to extend the 
present theory for circumstances of polydispersion by considering the recommended average 
droplet diameter. A detailed discussion of polydispersion effects is beyond the scope of the 
present investigation. Further measurements on both sound absorption and dispersion will be 
valuable in our understanding of the mechanisms and further development of the theoretical 
models for sound propagation in gas-vapor-droplet systems. 

5. Conclusion 

The effect of transpiration on sound propagation (attenuation and dispersion) in evaporating 
gas- vapor-droplet mixtures initially in equilibrium has been theoretically investigated. It is noted 



that at low droplet concentrations, two distinct absorption peaks due to mass transfer and 
viscosity-heat conduction effects are manifested, but at sufficiently high droplet concentration, 
significant interaction between the transport processes result in only one distinct attenuation 
peak. The results show that , transpiration has an appreciable effect on sound attenuation and 
dispersion at low and high droplet concentrations. For low droplet mass fractions (C^ »0.01) 
when mass transfer effects predominate over viscous and heat conduction effects, the strongest 
effect of transpiration on sound absorption is manifested at ryr, / = 1 , while the strongest 

effect of transpiration on sound dispersion is evident in the low frequency limit. At high drop 
concentrations, increases in the peak attenuation frequency and the dispersion coefficient are 
occasioned near « 1 . 
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Figure Captions 


Fig. 1. Comparison of dimensionless attenuation coefficient from Davidson's theory with the 
experimental data of Cole and Dobbins [14]. 

Fig. 2. (a) Comparison of dimensionless absorption coefficient from transpiration theory with 
data of Cole and Dobbins [14]. (b) Comparison of the dimensionless absorption coefficient with 
COT, from transpiration theory with the data of Cole and Dobbins [14]. 

Fig. 3. Predicted dispersion coefficient from transpiration theory for the test conditions of Cole 
and Dobbins [14]. 

Fig. 4 (a) Effect of transpiration on dimensionless absorption coefficient predicted by the present 
theory for high droplet concentrations, (b) Effect of transpiration on the dispersion coefficient 
predicted by the present theory for high droplet concentrations. 




Fig. 1. Comparison of dimensionless attenuation coefficient from Davidson's theory with the 
experimental data of Cole and Dobbins [14]. 




t m 


(a) 



0 ^ 


Fig. 2. (a) Comparison of dimensionless absorption coefficient from transpiration theory with 
data of Cole and Dobbins [14], (b) Comparison of dimensionless absorption coefficient with 
COT, from transpiration theory with the data of Cole and Dobbins [14] 
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Fig. 4 (a) Effect of transpiration on dimensionless absorption coefficient predicted by the present 
theory for high droplet concentrations, (b) Effect of transpiration on the dispersion coefficient 
predicted by the present theory for high droplet concentrations. 




